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Abstract 

Besides the standard quantum version of the Coulomb/Kepler problem, an 
alternative quantum model with not too dissimilar phenomenological (i.e., 
spectral and scattering) as well as mathematical (i.e., exact-solvability) prop- 
erties may be formulated and solved. Several aspects of this model are 
described. The paper is made self-contained by explaining the underlying 
innovative quantization strategy which assigns an entirely new role to sym- 
metries. 
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1 Introduction 



The list of the traditional roles of symmetries in Quantum Theory has re- 
cently been enriched by the increasingly active studies of the role of the so 
called PT— symmetries [1]. Typically, one considers a Hamiltonian composed 
of a kinetic- and potential-energy terms, 

H = -—A + V{x) 
2m 

and replaces the current requirement H = W of its Hermiticity in Hilbert 
space LF'iW^) := ']-i(f"^"'^''y^ by the less trivial physical Hermiticity (or rather 
"crypto-Hermiticity" ^) H = which must be constructively, ad hoc de- 
fined in anot/ier Hilbert space _ in the language of symmetries this 
means that, firstly, the "usual" Hermiticity H = W is reinterpreted as a 
formal time-reversal symmetry 

Hr=TH (1) 

(with a suitable time-reversal operator T). In the second step one rejects the 
constraint ([1]) as "too formal" ^ and postulates, instead, the Hermiticity of 
H = H* in an auxiliary Krein space with a suitable indefinite pseudometric 
V. In other words, in the language of symmetries one follows the recommen- 
dations of mathematicians |H |5] and re-facilitates the mathematics by the 
replacement of Eq. ([1]) by the modified requirement called VT symmetry, 

VrH = HVr. (2) 

In the context of physics the recipe finds its most ambitious theoretical en- 
couragement in the appeal of the concept of such a form of symmetry in 
relativistic quantum field theory [6], with V representing the parity and with 
T mimicking the time reversal in this implementation. 

Another explanation of the increasing popularity of the whole concept 
certainly lies in the amazing productivity of the PT— symmetrizations of 
various quantum systems. In the popularization of such a trick the key role 
has been played by the serendipitious letter [7] in which Bender and Boettcher 
proposed and demonstrated that the "PT— symmetrization of a given local 
Schrodinger equation in single dimension (i.e., a transition from constraint 
([1]) to constraint (|2]) at d = 1) may represent an efficient theoretical tool and 
way towards finding new phenomenologically useful Hamiltonians with real 
spectra. They also recommended to achieve this simply by the replacement 
of the original Hermitian potential V{x) by its non-Hermitian alternative 
V{ix) (thorough review offered by paper [3] may be recommended as an 
introductory reading). 

In our present paper we intend to return to the older application of this 
idea to the exactly solvable Coulomb /Kepler problem [8] and to its recent 
upgrades |[9j - [Hj. We intend to review the related recent theoretical devel- 
opments and to show how the incessant progress in the field applies to this 
particular but important example. 



2 The PT— symmetric version of the Coulomb 
problem 

In loc. cit., the standard quantum Coulomb/Kepler problem has been as- 
signed a new, non-equivalent quantum system version which is formally rep- 
resented by the PT— symmetric Schrodinger equation 

iZ 1 

+ — ^(x) = E^fx), L> — , 2L 4Z. 
X 2 

(3) 

We shall abbreviate here 2mE/h'^ = —k"^ and emphasize that even in the 
presence of the imaginary unit, this equation remains solvable in terms of 
the well known confluent hypergeometric functions, 

^,{x) = C^^^{x) + C2^2{x), (4) 

^^(x) = e-^''x^+\Fi{l + L + iZ/{2k), 2L + 2, 2kx) , (5) 

^2(x) = e~'=^a;-^iFi(-L + iZ/{2k), -2L, 2kx) . (6) 

The well known Coulomb - harmonic oscillator correspondence has been 
studied in [8] . In the role of a formal postulate it helped us to fix the physical 
asymptotic boundary conditions which would remain, otherwise, ambiguous 
(Ref. [12] and previous loci citati should be consulted for all details). 

It has been shown [13] that all of the models of the class sampled by 
Eqs. ([3]) + (H]) may be perceived, under certain conditions which we explained 
in Ref. [2], as fully compatible with the standard postulates of Quantum 
Mechanics. 

Our present paper will expose Eq. ([3]) as a special case of the broader class 
of Schrodinger equations which all exemplify an extension of quantum model- 
building strategies. In Ref. |2j we called this approach a "crypto-Hermitian" 
or "three-Hilbert-space" quantum mechanics. We shall emphasize here that 
the transition from the Hermitian to PT— symmetric language is extremely 
productive while, at the same time, its multistep nature often leads to con- 
ceptual misunderstandings. While "teaching by example", we shall try to 
clarify here some of the most blatant ones. 

3 The abstract formalism 

Several compact review papers [21 El El] may be recommended for reference. 
At the same time, an introductory explanation of the structure of Quan- 
tum Mechanics using the pseudometric in Krein space may be given, for our 
present purposes, a much shorter form. 

First of all the readers should be warned that in the majority of textbooks 
on Quantum Mechanics the meaning of the Dirac-ket symbols {i/j) E V is 
merely explained via examples. Most often one deals just with the most 
common quantum motion of a point particle inside a local potential well. 
Thus, it is assumed that there exists an operator Q of the particle position 
with eigenvalues g G M*^ and eigenkets \q). The standard (often called Dirac's) 
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dx^ 



x^ 



Hermitian conjugation is then represented by the transposition plus complex 
conjugation of any ket vector, yielding the bra vector, 

^{Dirac) . ^ _ (7) 

In the basis |g) one may then represent the ket of any pre-prepared state 
by the overlap = i.e., by the square-integrable wave function of 

the coordinate q. 

Our present purpose is not a criticism of this approach as such (interested 
readers may find such a criticism elsewhere [15]) but rather just of one of its 
consequences. In virtually all of the similar classes of examples, indeed, the 
vector space of states V is simply assumed endowed with the most common 
inner product 

/oo 
ra{q)i^b{q) dq (8) 
-oo 

with, possibly, the integration replaced by the infinite or finite summation. 
Thus, we may (and usually do) set V = i^^(M), etc. 

Without any real danger of misunderstanding we may speak here about 
the "friendly" Hilbert space of states "H^^-* = V, calling the variable q in 
Eq. ([8]) "the coordinate" . In parallel we usually perform a maximally conve- 
nient choice of the Hamiltonian H based on the so called principle of corre- 
spondence which encourages us to split the Hamiltonian into the kinetic and 
potential energies, H = T + V . Whenever the general interaction operator V 
is represented, say, by a kernel V{q, q') when acting upon the wave functions, 
this kernel is most often chosen as proportional to the Dirac's delta-function 
so that V becomes an elementary multiplicative operator V = Viocai = V{q)- 
Similarly, the most popular and preferred form of the "kinetic energy" T is 
a differential operator, say, T = Tiocai = —d'^/dq^ in single dimension and in 
the suitable units. 

The word of strong warning emerges when we perform a Fourier transfor- 
mation in "H^^^ so that the variable q becomes replaced by p (= momentum). 
One should rather denote the latter, Fourier-image space by the slightly dif- 
ferent symbol 1-L'^^\ therefore (with the superscript still abbreviating "phys- 
ical" [2]). 

Paradoxically, after the latter change of frame the kinetic operator Tiocai 
becomes multiplicative while Viocai becomes strongly non-local in momenta. 
Nevertheless, all this does not modify the overall paradigm. A truly deep 
change of the paradigm only comes with the models where the necessity of 
the observability of the coordinate q is abandoned completely. One may still 
start from the vector space of kets V but it makes sense to endow it with 
another Hilbert-space structure, via the inner product defined by an integral 
over a complex path, 

{^a\i^t) = J^ra{s)Ms)ds. (9) 

This is one of the most characteristic intermediate steps made in the so called 
PT— symmetric quantum theories [3]. The resulting loss of simplicity of the 
position operator Q changes the physics of course. The key point is that we 
lose the one-to-one correspondence between the integration path C and the 
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spectrum M of any coordinate-mimicking operator. The physics-independent 
optional variable s becomes purely formal. 

In such a setting our choice of the physical observables must still obey 
the old quantization paradigm, which is just set in a modified context. The 
loss of the observability of the coordinate proves essential, anyhow. For 
illustration one might recall the pedagogically motivated paper [16] in which, 
in a slightly provocative demonstration of the abstract nature of quantum 
theory, the variable s in Eq. (Q has been interpreted as an observable "time" 
of a hypothetical "quantum clock" system. 

Once we wish to understand our Coulombic Schrodinger eigenvalue prob- 
lem ([3]), we must make one more step and generalize further the inner prod- 
ucts ([9]). Such a second-step generalization of the inner product will certainly 
move us from the two Hilbert spaces H^^^ and Ti^^^ to the third one, viz., to 
the final and physical "standard" Hilbert space Ti^^^ (this notation is taken 
from Ref. 0). 

The introduction of the third Hilbert space forms the theoretical back- 
ground of an amendment of the traditional quantum mechanics, the key 
nonstandard features of which can be seen 

• in the admissibility of the complex potentials sampled by the power- 
law-anharmonic family V{s) = — (is)^"'"'^ of Ref. and generating the 
real and discrete bound-state spectra at any S > (cf. the proofs in 

my, 

• in the replacement of the usual real line of s by a complex curve 
C = C(R) which may even be, in principle, living on a complicated 
multisheeted Riemann surface [T2]; 

• in the theoretical imperative of the construction of certain operator 
(see below); 

• in the possibility of a systematic study of the discretizations and sim- 
plifications. 

In our present paper, the emphasis will be put on the last feature. 

At the stage of development where we did not yet explain the meaning 
and role of the operator (called Hilbert space metric) the theory remains 
incomplete. We already cannot rely upon a more or less safe guidance of 
quantization as offered by the principle of correspondence. Just a partial 
revitalization of such guidance is possible in the new context (cf., e.g., a nice 
example-based discussion of this point in Ref. [13]). 

This being said, the main theoretical obstacle lies in the vast ambiguity 
of the necessary appropriate generalization of the Hermitian conjugation as 
prescribed by Eq. ([7]). The general recipe (explained already in [IB] or, more 
explicitly, in [2]) is 9— dependent and reads 

^(general) ^ ((^|:=(^|e. (10) 

This means that using the language of wave functions ipls) with s G C we 
must replace the most common single-integral definition (jH]) of the inner prod- 
uct in the original "friendly" Hilbert space Ti^^^ by the more sophisticated 



double-integral formula 



(11) 
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In terms of an integral-operator- kernel representation B(s, s') of our abstract 
metric operator = ©"I" > this recipe defines the inner product which 
converts the same ket-vector space V into the amended and final, metric- 
dependent and physics-representing Hilbert space Ti^^^ of the very standard 
quantum theory (cf. [2J for more details). 

4 The upgraded formalism in applications 

In the attempted applications of all of the new ideas to the traditional bench- 
mark models like Coulomb scattering one may make use of its traditional 
merits (like, e.g., exact solvability) as well as of the flexibility of the choice 
of the PT— symmetric (i.e., complex and left-right symmetric) integration 
path, sampled in Ref. [10] as follows. 



This leads to new results of course. Typically, in spite of the non-unitarity of 
the scattering (remember that the Coulomb potential is strictly local!) the 
bound-state energies still emerge from the poles of the scattering matrix |10j . 

4.1 Discretizations 

The use of discretizations of the differential forms of Schrodinger operators 
may be, typically, Runge-Kutta-inspired. In practice, they are slowly becom- 
ing useful in sohd-state physics [I9], optics [20] and statistical physics [2Tj . 
Less expectedly, the use of lattice models proved crucial for the unitarity 
of the scattering. It has been shown [22] that the theory of scattering by 
non-Hermitian obstacles may be made unitary and consistent via a certain 
selfconsistently prepared transition to non-local potentials. Unfortunately, it 
is not yet clear how such a requirement of selfconsistency could be realized 
in the continuous limit of the discrete models. 

Whenever we discretize the coordinates and replace the differential Hamil- 
tonians by matrices with property H ^ H\ the above-reviewed theory ap- 
plies without changes. The usual Hilbert space becomes unphysical and it 
must be replaced by its unitarily non-equivalent correct alternative en- 
dowed with a sophisticated metric = 0*^"^^ 7^ / which defines the ad hoc 
inner product. 

In the discretized version of the theory, the integral kernel of the metric 
must merely be replaced by a matrix. Naturally, also the double integral f|TT]) 
gets replaced by the double sum. 




(12) 



N 



in 'H^^\ In this setting, the imaginary choice of the Coulomb couphng may 
still be made compatible with the standard postulates of Quantum Theory, 
provided only that it still generates the real, i.e., potentially observable spec- 
trum of the bound-state energies. 

Once we started our considerations from the imaginary Coulomb model 
defined along a continuous complex trajectory, we may expect that many 
of its properties will survive also the transition to its discrete descendants. 
For inspiration we may recall Ref. [10] where the bound-state energies were 
shown to coincide with the poles of transmission coefficients. Still, as long 
as the potential 7^ VMs local, the unitarity of the scattering cannot be 
required (TUl 122] • At the same time, the unitarity of the time evolution of 
the system itself may be achieved. Indeed, although the Hamiltonian H is 
non-Hermitian in 'H^^\ (abbreviated H 7^ H^)^ it is Hermitian in 'H^^\ This 
feature is called cryptohermiticity, requiring H = alias 

Here, the operator or matrix Q^^^ is precisely the one which defines the 
physical inner product. 

This being said, the loss of easy constructions is a problem [23]. Still, the 
discretization of the coordinates may be recommended as the recipe. 

4.2 Interpretations 

The generic VT— symmetric quantum model describes a closed system de- 
fined via a doublet of operators H, 6 or via a triplet of operators (adding a 
new observable A and having, typically, Hamiltonian H 7^ accompanied 
by a charge [3]), etc. In other words, the dynamical content of phenomeno- 
logical quantum models is encoded in Hamiltonian H and in metric G. In 
this setting the metric G guarantees the unitarity of time evolution in an 
ad hoc, "standard" Hilbert space [18], to be denoted by the symbol Ti^^^ in 
what follows. In addition, one can also impose some other, phenomenologi- 
cally motivated requirements like a short-range smearing of coordinates [23], 
etc. 

One of the remarkable features of such an upgrade of applications of 
quantum mechanics may be seen in the robust nature of its "first principles" 
which remain unchanged. Thus, its traditional probabilistic interpretation 
is not changed (notice that it practically did not change during the last cca 
eighty years!). In the language of textbooks one could speak just about the 
use of a non-unitary generalizations ~ \/G of the most common Fourier 
transformations. Still, the new physics behind the trick may be nontrivial 
(in nuclear physics, for example, the mapping fl (called Dyson's [18]) was 
used to represent fermions as images of bosons) . 

Among the most innovative consequences of the upgraded formulation 
of quantum models one notices, first of all, the existence and possibility of 
constructions of a horizon [2S]. Formally, this notion coincides with the set 
&D of the Kato's [26] exceptional points in the (real or complex) manifold 
V of available free parameters (like coupling strengths, etc). The practical 
appeal of this notion may be based, e.g., on its tunability [2Z] and/or a new 
physics near instabilities and quantum catastrophes [2B]- 



As another emergent concept one should hst fundamental length, i.e., 
a quantity 6 defined, in the simplified discrete models, as the number of 
diagonals in the metric which is required to possess a band-matrix form, 
6mn = for \m — n\ > 6. In this context one might mention the first papers 
devoted to the study of PT— symmetric quantum graphs [2S] in which one 
might search for a connection between the fragile parts of the spectrum and 
the topological characteristics of the underlying graph structure. 

Last but not least, it is necessary to emphasize the challenging character 
of a generic scenario with more observables, each of which may be responsible 
for its own part of the physical horizon, "invisible" from the point of view 
of the other observables. In other words, a lot of work is still to be done 
before one could speak about a "classification" of exceptional points (i.e., 
about a a sort of "quantum theory of catastrophes") - the first attempts 
in this direction only dealt with the hardly realistic, too oversimplified and 
schematic quantum systems [5U] . 

5 Coulomb potential V{xj) = i/xj 

The main weak point of the above-cited choice of the Coulomb potential 
may be identified not only with its strict locality (i.e., with the necessary 
loss of the unitarity of the scattering, cf. also the detailed study [23] in this 
respect) but also with the difficulties encountered during transition to any 
model which would not be exactly solvable. For this reason, our present use 
of the Runge-Kutta-inspired discretization will help also in the case of the 
Coulomb potential. 

As long as our present main ambition is the presentation of the upgraded 
formalism, we shall try to simplify many inessential mathematical aspects of 
our Coulomb/Kepler model. In parallel, we shall also try to treat this poten- 
tial as a special case of a broader class of forces. For the sake of definitness 
and in a way insspired by Ref. [?], we shall pick up the class V{x) = ix^ 
with a real exponent which does not lie too far from its Coulombic value of 
z = -1. 

For this purpose we must replace, first of all, the typical differential 
Schrodinger equation 

- — Mx) + V{x) 4j{x) = En ^n(x) ^(±A) = , 1< A < oo (13) 
ax'^ 

by its discrete version (i.e., approximation or analogue) 

An equidistant grid of the Runge-Kutta points Xk = —A + k h with k = 
0, 1, . . . , -|- 1 will be used. In this sense, also the standard general double- 
integral inner product will be replaced by the above-mentioned double sum, 
etc. 

Naturally, the discretization recipe also involves the change of the asymp- 
totic boundary conditions, with x^+i = A, h = 2A/{N + 1) and V'(%) = 
ip{x]\f+i) = 0. In other words, the eigenvalue calculations become reduced to 
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the mere diagonalizations of the by matrix Hamiltonians 



( 2 + /i2\/(xi; 
-1 



V 



-1 

2 + h^V{x2) 
-1 



-1 

2 + h^V{x^) 



-1 2 + /iV(x^) / 



(15) 

The insertion of any potential V{xj) = ixj will lead to the eigenvalues ej : = 

h'^Ej^^ G (0, 4) which must be computed numerically in general. 

In some applications the transition to the continuous limit = oo is 
made or, at worst, postponed till the end of the calculations. In the present 
methodical context we shall rather keep the dimension constant and, in 
fact, not too large. 

The insertion of formula V{xj) = i/xj in Eq. (fT5|) with even A^ = 2K 
yields the sequence of the discrete PT— symmetric Coulomb Hamiltonians 



2-ia/{2K -1] 
-1 








2-ia/3 



-1 



-1 



2 + ia 



-1 



-1 2 + ia/3 



In its first nontrivial example let us set A^ = 4, 



2 + ia/{2K-l) 
(16) 



ia3^ 
-1 





-1 
2 — ia 
-1 






-1 

2 + ia 
-1 







-1 

2 + 2 a 3" 



-1 



(17) 



We see that a natural generalization may be targeted not only at the growing 
dimensions A^ > 4 but also towards the small deviations of the exponent z 
from its Coulombic value. Empirically, one can verify that in both of these 
directions, the spectral loci (i.e., eigenvalues e^^'^\a)) remain topologically 
the same. More precisely, at a fixed A^, the topology of the Coulomb-potential 
pattern as sampled by Figs. [T]-[3]may be expected to survive all the negative 
exponents z 

At N = 4 the model is exactly solvable at the Coulombic exponent z = 
— 1. The secular equation 



E^ + 



10 . 
21 + — 
9 



E^ + 



40 



20 E + 5 + l/9a^ + 5a^ 







n 



generates the closed-form spectrum 



e{a) = 2 ± 1/6 V 54 - 20 a2 ± 2 V405 - 720 + 64 



which is real iff \a\ < 3/4 V 10 - 4 ^ 0.7706147226. 

At = 6 the model is still exactly solvable at z 
secular equation 



- 12 E^ 



55 + 



259 
225 



E* + -120 



2072 
225 



■1, yielding the 



E^+ 



5894 2 7 4 > 2 

+ 126 + + — aM ^ + 

225 45 



-56 - 



280 



28 
45 



E+ 



+7+Ua^ + -a^ + — 
9 225 







which may be solved using Cardano formulae. Although the closed form of 
the spectrum becomes extremely clumsy in this representation, it decisively 
facilitates the graphical representation of the spectral loci which all appear 
topologically equivalent to the vertical array of circles. In particular, the 
survival of the exact solvability of the problem enables us to conclude that 
the whole A^ = 6 spectrum remains real iff |a| < 0.589586. 

These results indicate that the topological pattern remains generic and 
A^— independent. Such a conjecture is persuasively confirmed by the larger— A^ 
graphical samples which are presented in Figs. [2] and [31 




Figure 1: The a— dependence of energies e{a) aX N = 4. 



5.1 Graphical methods 

Numerical evaluation of the spectra is sampled in Figs. [H |2] and |3l These 
graphical constructions indicate that at any A^ = 2K the spectrum is real iff 
a e and fully complex for a i . where a^^^) 

is a quickly decreasing function of K. 

The latter observation may be interpreted in two ways. For the finite 
lattices in which the numerical value of parameter a is fixed, the reality of the 
spectrum is, undoubtedly, fragile. In the alternative approach in which our 
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model serves just as a simulation of the (NB: exactly solvable!) differential- 
equation system, the definition of parameter a is prescribed by the Runge- 
Kutta recipe (see above). For this reason, its numerical value decreases, with 
N ^ much more quickly than a*^^-*. This implies that in the latter setting the 
reality of the spectrum may be declared robust and guaranteed. 




-1 -0.5 0.5 1 



Figure 2: The a— dependence of energies e{a) at = 8. 

Marginally, we may add that in the former scenario using small and fixed 
A^, the loss of the reality of the spectrum is caused by the confluence of the 
ground state with the first excited state and by their subsequent complexi- 
fication. Due to the up-down symmetry of the spectrum, this instability is 
paralleled by the upper two states of course. 
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Figure 3: The a— dependence of energies e{a) dX N = 14. 



6 Hermit izat ions 

6.1 The N = 2 metrics 

The condition of hidden Hermiticity of any Hamiltonian with real spectrum 
is often called Dieudonne equation [2], 

HiN)^ [Hi^)f:=Q-^[H^^)]^Q. (18) 



Its structure is best illustrated at iV = 2 where the metrics form just the 
two-parametric family 



(k,m) 



k km — ika 
km + ika k 



k,m & 



(19) 

With the condition a G (—1, 1) of the reality of energies and with the metric- 
positivity constraint 



e = e± = k± VkW + k^a^ > 



(20) 



we may conclude that k must be positive and larger than the square root. We 
may reparametrize a = cos/3 sin 7, m = cos/3 cos 7, /3 G (0,7r) and 7 G (0, vr) 
and — 1 < cos /3 < 1 and get the final result 



e = e{i/(2)[a(/3,7)]} 



[fc,m(/3,7)] 



k- 



1 e 

e^'^ cos /3 



-'T cos /3 
1 



(21) 



In a search for the ot/ier eligible observables with crypto-Hermiticity property 

^(/3,7) 0[fc,»n(/3,7)] = 0[fc,m(,3,7)] A(/3,7) • (22) 

the use of the ansatz 

^ G + ig B + ib 
C + ic D + id 



A 



(23) 



leads to the four real constraints imposed upon eight free parameters. The 
family of observables is four-parametric, therefore. Three constraints define 
B, C and G — D. The remaining one relates the sums cs = b + c and 
9y, = 9 + d and leads to the unique solution (^s = 0. We may conclude that 
from the input m = m{/3, 7) and a = a{f3, 7) one gets the class of admissible 
observables 



A = AiD,b,c,g) = - 



Da — b — c + iga, g — bm + iba 
g + cm + ica, Da — iga 



(24) 



In particular, the initial Hamiltonian is reobtained at D = 2, G = 2, b = c = 
0, B = C = —1 and g = —a (= —d). 

In the literature the concept of charge is considered particularly useful [3]. 
Its essence lies, in the present model, in an additional auxiliary assumption 



[H]^V = VH 



(25) 



where V is the operator of parity. Under this assumption a unique metric is 
sought such that a very specific metric called CVT metric is prescribed by 
formula G)('^^'^) = CV where C is called "charge". 

Thus, at a given N we may define the parity V = V^^^ which con- 
tains units along the secondary diagonal, i.e., Vm,n = liffm + n = A^ + l 
while Vm)i = otherwise. Incidentally one may note that p(^) = limit of 



2)( 



such that ^ 0, m^^) ^ 00, A;(^)m(^) ^ 1. 



The key merit of the use of charge is that its use makes the metric unique, 



c(^^'^) = k 



-ia 1 
1 ia 



(26) 



Moreover, it also represents one of the special cases of observable A using 
D = b = c = and g = — cos/3/ sin/3 = —y/k'^ — 1. Indeed, from Q^'^'^'^ = 
CV we have 

u V 



c = e[ffp>( 



y z 



(27) 



yielding v = y = k and z = u* = ke^'^ cos /3. Then, condition = I requires 
that 7 = 7^^^^) = 7r/2 (i.e., a = cos/3). Thus, we have /3 = /3('^^^) such that 
sin /3*^'^'^^) = l/k. This proves the above statement. 



6.2 The N = 4 metrics 

With the natural ansatz for G ^H^^'> (a, z 



k m — ikw 

m + ikw r 

W h + i {kw + ra) 
Z W 



{k,m,r,h) 

w* 

h — i {kw + ra) 
r 

m + ikw 



Z* 
W* 
m — ikw 
k 



{21 



where w = w{z, a) = 3^a and 

W = W{k, m, r) = —w'^k + r — k — kwa + i {wm + ma) , 

Z = Z{k, m, r, h) = mc? — w'^m — m + h — i {kw — ka — kwa^ — rw + w'^k) , 

the problem of the determination of the domain of positivity of the metric 
starts to be merely tractable graphically. 

For the numerous practical purposes the metric is sought in a special 
form. One of the phenomenologically inspired options is the choice of the 
matrix with units along its main diagonal, k = r = 1. In addition, let us 
select m = h = and compute W{1,0,1) = —w{w + a), Z(l, 0,1,0) = 
i{a + wa^ — w^). Such a restricet construction leads to the following four 
closed- form eigenvalues of the metric [if ( 



«' I (1,0,1,0) '^1^' 



= 1 + ^ {w-a^w + w'^)±]^ 



e^ = l- \ {w - a^w + w^) ±]- 



A 



± 



w 



+ (2 - 2 a^) ti;^ + (±8 + 4 a) + (5 ± 8 a + 6 + a*) w'^+ 



+ {Aa + Aa^) w + Aa^ . 



Some details and numerical results of its analysis may be found elsewhere 
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7 Discussion 



In a climax of our present discussion of the discrete PT— symmetric Coulomb 
problem characterized by the purely imaginary coupling constant, let us now 
summarize the overall method via the following scheme 



textbook level quantum theory : 
prohibitively complicated Hamiltonian f) 
generating unitary time evolution 

trivial 



P: physics 



calculations = practically impossible 



simplification 



unitary equivalence 



state i\i is represented 
in the false Hilbert space 



I F: I calculations = feasible 
physical meaning = lost 
H = non — Hermitian 



hermitization 



amended inner product 



S: 



standardized representation 
picture = synthesis 
physics = reinstalled 
H = Hermitian 



which characterizes the "three-Hilbert-space" pattern of quantization as de- 
scribed in Ref. |2j as a recipe in which the usual Schrodinger equation 

H\ilJn) = E^\iljn) (29) 

finds the standard probabilistic interpretation even if the Hamiltonian matrix 
(with real spectrum) proves manifestly non-Hermitian. 

The specific feature of non-Hermitian matrices H may be seen in their 
ability of having the reality of their spectra controlled by a parameter (for this 
purpose we used a in our present models). In other words, one can simulate 
the abrupt loss of the stability of the time evolution of the system by a mere 
smooth change of this parameter. In other words, we may speak about a 
non-empty (quasi-) Hermiticity domain of parameters, with the qualitative 
changes of physics at its boundary, and with a guaranteed reality of the 
spectrum in its interior. 

In our present paper we emphasized that another important aspect of 
physics with real spectra but non-Hermitian matrices of observables lies in 
the necessity of a fine-tuning mediated by the Dieudonne equation. 

Typically, a given Hamiltonian H must be assigned a Hermitizing metric 
9. As long as we merely considered < oo, we could avoid any difficulties 
by simply solving the second, conjugate Schrodinger equation 

{{^|J^\H = E^{{^Pm\ (30) 
which may be also written in the form 

This enabled us to work with the solutions as forming a bicomplete and 
biorthogonal basis, 

N-l ^ 

^ =Y1 1^"^ //I I / \ ^^^'^l ' ((^mlV^n) = S^,n {{i^nli^n) ■ (31) 

„=0 \Wn\Wn) 



The main benefit may be then found in the closed formula 



N-l 

|2 



= E 1^-)) l'^ 



n \\ Yn 



n=0 



which defines all of the eligible metrics. This, in its turn, specifies all the 
dynamics given by the operator doublet (if(A), 6(/t)). 
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